Phase estimation with photon number constraint 
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Many researches proposed the use of the noon state as the input state for phase estimation, which 
is one topic of quantum metrology. This is because the input noon state provides the maximum 
Fisher information at the specific point. However, the Fisher information does not necessarily give 
the attainable bound for estimation error. In this paper, we adopt the local asymptotic mini-max 
criterion as well as the mini-max criterion, and show that the maximum Fisher information does 
not give the attainable bound for estimation error under these criteria in the phase estimation. We 
also propose the optimal input state under the constraints for photon number of the input state 
instead of the noon state. 



I. INTRODUCTION 

As one of the most simple methods in quantum 
metrology .phase estimation has been treated by many 
authors [3-0| and its experimental demonstration has 
been reported by many groups [8HTlj|. As is pointed by 
[2lj . the phase estimation with Mach-Zehnder (MZ) in- 
terferometer is important from the view point of the stan- 
dard configuration. Many experiment groups realized the 
noon state. The reason is why the noon state provides the 
optimal Fisher information. In the independent and iden- 
tical distributed setting of state estimation, Fisher infor- 
mation provides the asymptotic bound of mean square 
error (MSE) [l7l - fl9{ . However, the phase estimation does 
not fall in this setting. Hence, it has not been discussed 
sufficiently whether the minimum estimation error is re- 
alized when the noon state is inputed. 

In this paper, we treat the phase estimation with two 
kinds of photon number constraints. One is the con- 
straint for the average of the square of the photon number 
concerning the input state. The other is the constraint 
for the maximum photon number concerning the input 
state. In order to treat the phase estimation, we focus 
on two criteria. One is the minimum mean square er- 
ror (MSE) under the locally unbiased (LUB) condition, 
which coincides with the SLD Fisher information 

EMI- 

The other is the mini-max value of the MSE, which mini- 
mizes the worst MSE. The latter is essentially equivalent 
with the optimization with the covariant restriction for 
measurement [HI . As an intermediate concept, we focus 
on the local asymptotic mini-max criterion, which was 
introduced by Hajek [l2| and was applied to quantum 
channel estimation by Hayashi Then, we can define 
three kinds of bounds under two kinds of constraint. The 
main purpose of this paper is the comparison between 
these three bounds with both constraints. The second 
purpose is to check the validity of use of the noon state 
as the input state and seeking the optimal input state. 

The organization of this paper is the following. In 
section UIl we treat the phase estimation problem under 
the locally unbiased condition. In section IIII1 we clar- 
ify the relation between the single-application case and 
the multiple-application case. In section HVl we discuss 



the locally asymptotic mini-max criterion. In section fVl 
we treat the global mini-max criterion and the relation 
among three kinds of criteria under the asymptotic limit. 
In section [VTI we treat the MSE of the noon state with 
both of the global and local asymptotic mini-max crite- 
ria. 



II. SLD FISHER INFORMATION AND 
LOCALLY UNBIASED CONDITION 

In the two- mode photonic system L 2 (]R 2 ), the phase 
operator Ug is given by 



1 71, m)(n, m\ 



n,m>0 



with 9 6 [— 7r, 7r] . When the true parameter 9 is unknown, 
we can estimate the parameter by inputing the known 
state \<j>) and measuring the output state Ug\4>). Hence, 
when the measurement corresponds to the POVM M = 
{Mg}, our estimator is described by the pair of the input 
state |</>) and the POVM M as Fig[TJ 

Input Unknown phase Output Measurement estimate 
state operator state 

k) ' 1^1— ~ U eW) ~~\M\ 



FIG. 1: Scheme of phase estimation 

In this case, the estimate 9 obeys the distribution 
(cj)\UgMgUg\(j)). The mean square error (MSE) is given 
by 

E e (\4>),M):= f (9 - Te(6)) 2 (<t>\UtM § U e \<f>)d§, 

J — 7T 

where 

( 9 + 2tt ii9 <9-n 
T 9 {9) := I § tf9-ir<9<9 + ir 

{ 9 - 2tt if 9 + 7T < 9. 
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When the input state U$\<p) is given, our problem is es- 
timation of the parameter 6 with the pure state family 
{U$\<p)((/)\Ug}. In order to treat the lower bound of the 
MSE, we focus on the symmetric logarithmic derivative 
(SLD) Fisher information (20j 



4|»0#| 2 . (1) 



This equation can be shown by the following way: 
^Ug\<f>){<f>\Ul = X|0><0| + |0)<0|Xt, where X = 

% ~ (<t>\^\<f>)- Thus, jgUeimWt = 2(AT|0)<0| + 
|0>(0|Xt) o |0) where A o B = AB+ba Theilj 
J e (|0)) = Tr(2(V|0) (0| + |0) (0|Xt )) 2 10) (0| , which equals 
the right hand side of ([T]). 

In statistical inference, an unbiased estimator plays an 
important role. However, since the parameter space is 
the interval [— 7r,7r], it is difficult to define the unbiased 
condition. It is possible to define the locally unbiased 
condition. When the estimator (|0), M) satisfies the con- 
dition 



T e {e){<t>\UlMfJe \<t>)de = 6 Q 



T e {ff) — {(j>\UlM 6 Ue\<i>)\8=eodB = 1, 



it is called locally unbiased at 9q. When an estima- 
tor (|0), M) is locally unbiased (LUB) at , the SLD 
Cramer-Rao inequality 



E e (\4>),M)> J e (|0))- 



(2) 



Its proof is similar to the usual case of the SLD Cramer- 
Rao ine quality [3-G1] because the Schwarz inequality 
plays the same rule. 

In this problem, when we input the state with a larger 
photon number, the MSE becomes smaller. So, it is 
suitable to constraint the photon number of the input 
state. As the first kind of constraint, we restrict the 
average of the square photon number (N 2 ) av (\(j))) :— 
(0|iV 2 |0) of the input state |0), where N := J2 n m( n + 
m)|n, m)(n, m\. As the second kind of constraint, we re- 
strict the maximum of the photon number N max {\<p)) := 
max„ jTO {n + m\{n,m\4>) ^ 0} of the input state |0). So, 
we define the two bounds of MSE as follows: 



G a v,e{E) 
:=inf 1^(10), M) 

C ma xfi{E) 

:=inf (i0(|0),M) 



(|0),M)is LUB at O 

(iv 2 UI0))<£ 2 



(|0), M)is LUB at Q 
N max (\cj))) < E 



Since J e (|0)) < 4(0|iV 2 | 
equality ([2]) implies 



the SLD Cramer-Rao in- 



C r naxfi{E) > C aV: g(E) > . 



When E is an integer n and the input state |0) is the noon 
state |0 n , noon ) := ^flra.O) + |0,n)), Jed&n^oon))' 1 = 
j^2. There exists a POVM satisfying the locally unbi- 
ased condition whose MSE attains the inverse of the SLD 
Fisher information. Hence, 



Cmax,e(n) = Cav,e(n) 



1 

4n 2 ' 



In the state estimation, it is known that there exists 
a sequence of estimators that asymptotically attains the 
inverse of the SLD Fisher information at all points (TtI— 
[l9j . Hence, we expect that the bounds 



a 



max, e nm E Cmax,e(E) — - (3) 

E— ¥00 4 



c„ 



:= lim E 2 C avfi {E) = - 

E— >oo 4 



(4) 



can be attained asymptotically in all points by a single 
sequence of estimator under the respective restrictions. 
However, the same argument does not necessarily hold 
in the case of channel estimation. 



III. MULTIPLE APPLICATION 

As an extension of the formulation given in Section 
ITTl we assume that the m-fold tensor product system of 
L 2 (K 2 ), i.e., L 2 (R 2 )® m , is allowed as the input system 
and the unknown operator Uf m is applied. It is required 
to estimate the unknown parameter 9. In this case, we 
consider the number operator N m := J^kLi I® k ~ 1 ® N <£> 
j®m-k ^ anc j estimator is given by the pair of the input 
state |0) on L 2 (R 2 )® m and a POVM M on L 2 {M?)® m . 
Now, we define the projection 



P a .b ■= ^2^2\ai,h){ai,h\ 

#1 #2 



■ <8 \a r . 



,b m }(a 



tni b rn I : 



where #1 is a\, . . . , a m : J2j a j = a an d #2 is &i , . . . , b m : 
bj = b. Then, we define the normalized vec- 
tor \<t> atb ) := pr-^yj Pa,b\4>) and the isometry V m ■= 
E<M>M>(<?W>I from L 2 (R 2 )® m to L 2 (R 2 ). Then, the 
pair of the state V m \4>) and the POVM {V m M§V^} gives 
an estimator on the single input system L 2 (R 2 ). The 
outcome of this estimator has the same statistical be- 
havior of the original estimator (|0),M) on L 2 (R 2 )® m . 
Further, the average of the square number of V m \4>) 
equals that of |0). Conversely, any estimator (|0), {Ms}) 
on the single input system L 2 (R 2 ) can be simulated by 
an estimator (W^|0), {W^MsV m }) on L 2 (M 2 )® m , where 
W m := ^ a b \ip a ,b)(a, b\ and \4> a ,b) is a normalized vector 
on the range of P a ,b- So, the optimization of estimator 
(|0), {Ms}) on the single input system L 2 (R 2 ) with con- 
straint concerning average of square of photon number is 
mathematically equivalent with that on the m-fold ten- 
sor product system L 2 (R 2 )® m with constraint concerning 
average of square of photon number. 
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Next, we consider the case when only state in the sub- 
space spanned by {|00), |10), |01)} is allowed as the input 
system. In this case, the input system is 3-dimensional 
and is written by %. Its m-fold input system is given 
as H® m . For a given input state \<j>) € H® m , the state 
V m \(j)} satisfies the condition max 0i i,{o + b\(a, b\V m \(ft) ^ 
0} < m. Conversely, when the input state \4>) satis- 
fies the condition max aj j,{a + b\{a, b\(j>) 7^ 0} < to, and 
when \ipa,b) is chosen from the intersection the range of 
P a , b and H® 171 , the state W m \cj>) belongs to U® m . The 
estimator (W m \(f>, {WmMgW^}) simulates the estimator 
{Mg}). So, the optimization of estimator on the 
single input system L 2 (R 2 ) with constraint concerning 
maximum of photon number is mathematically equiva- 
lent with that on the m-fold tensor product system %. 

Overall, multiple application setting is a physical sit- 
uation different from single application setting. This is 
because multiple application requires multiple applica- 
tion on the same system L 2 (R 2 ). However, this setting 
can be mathematically simulated by single application 
setting. This kind of simulation plays an important role 
in the next section. 



IV. LOCALLY ASYMPTOTIC MINI-MAX 
CRITERION 

In this section, we consider the locally asymptotic 
mini-max criterion. In statistics, we often treat mini- 
max criterion, in which, we optimize the worst case. In 
the following, we optimize the worst MSE among the e- 
neighborhood of the given point 9$. We treat the asymp- 
totic behavior of this optimum value and consider the 
limit e — > concerning the asymptotic coefficient. That 
is, we define 



Cavfi 

C a v,e{E, e) 
Ce {\4>) :e) 

Cmax : 9 

9(E,e) 



Aim lim E 2 C aVi g(E,e) 

e-)-0_E-i-oo 

-M{C e m:e)\(N*) av m<E 2 } 

m 

:inf sup E e (\4>),M) 
M eeu(e,e Q ) 

Aim lim E 2 C max g(E,e) 

£->0£;->oo 

= tof{C fl (|0) ■t)\N max (\4>))<E}. 

m 



We also define similar values with multiple application. 

C a vfi Q ,rn 

:= inf { sup Eg{\cj>),M)\{Nl l ) av {\ ( l>))<E- 2 m 2 } 

i\<t>),M) o e u(e,9 ) 
Cmaxfio^m 

:= inf { sup E e (\4>),M)\ \(f>) G H® m }, 

(\4>),M) 6<EU(efi ) 

where (|</>),M) is an estimator on L 2 (R 2 )® m . As is dis- 
cussed in Section UTTl estimators on L 2 (M 2 )® m can be 
simulated by estimators on L 2 (R 2 ) with large photon 



number states. Hence, 

C„„j =lim lim m 2 E 2 C a vfi,m(l, e) 

e— ^0 m— foo 

Cmax, e = hm lim m 2 C ma x,e ,m< 

e— >U m— >oo 

Further, Proposition 1 of [l| implies that 



lim lim to C av a m (l, e) < C av .e 

e— m->oo 

lim lim m C max a m < 



>0 m—toc 



That is, 



Cavfi ^ C av fi 
Crnax.O — Cmax. 



V. GLOBAL MINI-MAX CRITERION 

Next, we consider global optimization, i.e., global mini- 
max criterion. For this purpose, we consider group co- 
variant condition for our POVM M. In this optimization 
with photon number constraint, we can restrict the input 
state on the subspace Hi spanned by {|n,0)} U {|0, n)}. 
A POVM M — {Mg} is covariant when there exists a 

vector |T) such that M § d§ = ±U § \T)(T\Utd9, where 



\T) = El 



|0,n), and 9 n and 9 n 



are arbitrary coefficients. As is shown by Holevo [13j, 
since the stochastic behavior of the error 9 — 9 of the co- 
variant POVM does not depend on the true parameter 
9, any input state |</>) satisfies 

inf max E 9 (\6),M)= inf max Eg(\6),M), 

M 8 M-.COV 9 

Hence, when we optimize the worst value 
maxa Eg(\<p), M), it is sufficient to treat the covari- 
ant POVM. That is, our estimator is written by the pair 
of input state |</>) = Y^=o a n\n,0) + J2^=i a -n|0, n) 
and |T) = £~„e*»|»,0) + £~=i e ifl »|0,n>, where 
a = {onj^-oo G L 2 (Z). This estimator is 
simulated by the pair of the input state \(f>') = 



Eoo 
71=0 e 



-i0, 



a n n 



>0) + £„=i e ™a_„|0, n) and the co- 



variant POVM given by \T ) = J2Zo K 0)+E„=i |0, n). 
This is because \{T\u\Ue\(l))\ 2 = \{T \u\UeW)\ 2 ■ Hence, 

u u 

it is sufficient to treat the case with |T) = \Tq). In 
this case, we define a' n := e~ l9n a n for n > and 
a'_ n := e~ n a_„ for n > 0. Then, when the error 9 — 9 
obeys the square \F~ 1 {a'){9)\ 2 of the inverse Fourier 



transform F^ 1 (a!)(9) 
series of a' = {a' n }^ = _ { 
Thus, we obtain 



Eoo 
n— — c 



-in9 



of discrete 



C{\4>)) := inf sup Eg(\<p),M) 

M a 



inf { — 

o':|a'|=|a n | 27T 



ij-vx^i 2 ^ 
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and 



C av (E) 

-- inf {sup E e (\<P), M)\(ct>\N 2 \cP) < E 2 } 
(\<PM) e 



■■ inf { — 

a:||a|| = l 2lT 



Cmax \E) 



9 2 \^-\a)(9)\ 2 de\(N\ v (a)<E 2 } 

(5) 



: inf {sup Eem^imN 2 ^) < E 2 } 
<,\4>M) e 



■■ inf { — 

a:||a||=l 2n 



\T-\a)(6)\ 2 d6\N max (a) < E}, 



(0) 



(7) 
(8) 
(9) 



where (N 2 ) av (a) := J2 n n \ a n\ . N max (a) 
max„{|n| \a n ^ 0}, and ||a|| 2 := Er=-oo l a «l 2 - 
From the definitions, we obtain 

cm > com ■■ c) 

C av (E) > C av ,9(E,e) 

C m ax(E) > C maXi g(E,e), 

and the asymptotic limits are defined by 



C av := lim E 2 C av {E) 

E— >oo 

Cmax hm E 2 C max (E). 

E^-co 



Given an arbitrary sequence = {a_E,n}^L_oo of ele- 
ments of £ 2 (Z), we focus on the function 

f(x) := lim \fEa E ,E x . (10) 
Then, the relations ([5]) and ([5]) imply that 



Cav = f Mf\ p2 \f)\(f\Q 2 \f)<l} (11) 



inf 

/eL 2 ([-l,l]):||/ll=l 



(/l^l/). 



(12) 



where Q and P are the position and momentum operators 
on L 2 (R). The equation ([T2"|) is the same as the result 
obtained by Imai and Hayashi Q. 

Next, we focus on the specific sequence of the input 
states \4>e) whose coefficient satisfies (fTO)) . Similar to 
(|TT|) and ((T2l), we obtain 



lim E n C{\4> E )) 

inf (g\P 2 \g). 

The uncertainly relation implies that 
C - 



(13) 



and the infimum is attained when f(x) = e t. The 
latter is given by 

and the infimum is attained when f(x) — 
(27T) 1 / 4 sin ElidgJ t Comparing the right hand sides 
of © and © with those of (IT]) and jU]), the bounds 
C av and C m ax can be attained by sequence of input 
states 



ci(E e ^|n, 

n=0 



o> + E 



e 



|0,n» (14) 



n=l 



and 



A . TrQE + n + 1) , 
c 2 (2_^ sm — ; — |n, 0) 



71=0 



2E + 2 



E 



. 7r(£ + n+l) 

3m 

2E + 2 



|0,n», 



(15) 



where c\ and c 2 are normalized constants. 
From the relations ((8|) and (O, we obtain 

Cav ^ Cav, 9 
Cmax -^L CmaxM- 



(16) 
(17) 



As is shown in Proposition 3 of Hayashi [1], the 
bounds Cavfi = ^e^o^ m ^oom 2 E 2 C a vM, m (^,E) and 
Cmax,e = lim ( ^ lim^oo m 2 C mal ,« „ can be globally 
attained by two-step methods [l?], [l8[ asymptotically. 
Since any estimator on the multiple system L 2 (R 2 )® m 
can be simulated by a suitable estimator on the single 
input system L 2 (R 2 ), the opposite inequalities of (fT6|) 
and (fl~7| hold. That is, combining (j3|) and ((4|), we obtain 



Cav, 9 — C a v — ^ — C a v,9 



c ~?-> 1 



(18) 
(19) 



So, when we adopt the maximum photon number con- 
straint, there is a non- negligible difference between the 
local asymptotic mini-max bound and the bound with 
locally unbiased condition. 



VI. ANALYSIS ON NOON STATE 

Now, we consider the asymptotic limit of the perfor- 
mance with the locally asymptotic min-max criterion 
when the noon state -^(|n,0) + |0, n)) is inputed and 
the covariant measurement is applied. In this case, first, 
we fix an arbitrary small real number e > and focus 
on the neighborhood C/(e,#o) for an arbitrary point 9o. 
When n > 2n/e, the relation U g ^(\n,0) + |0,n)) = 
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Ue + 2ir/n^(\n,Q) + |0, n)) holds and 9 and + ^ belongs 
to U(e,6o). That is, we cannot distinguish two parame- 
ters 9 and 9 + ^ . Even if we could estimate the unknown 
parameter 9 mod — perfectly, it is not easy to estimate 
the parameter 9 in the parameter space f/(e, #q). In the 
following, we consider the case of 9 — 9 — e mod — . As- 
sume that we could the unknown parameter 9 mod — . 
In this case, there are still K := [^rj candidates with the 
width — in the parameter space U(e,9 ). Assume that, 
we decide the true parameter to be ^-i + 9o — e with the 
probability pi . Then, the mini-max error is evaluated as 
follows. 

K 



1=1 



K 



2tt " K 

= (— ) 2 max (j 2 - 2j V i Pl + V i 2 Pi ) 
n j=i k L — ' L — ' 

i=l t=l 
27T ^ K K 

=( ^ )2( j = I ? ax K ( - 7 ~ S *^ )2 + Z * 2k ~ E zpi)2) 

i— 1 z— 1 z— 1 

2tt ^ 

>(^) 2 (^/2) 2 = (^L-J) 2 . 

n n 7r 

That is, we obtain 

C, o (i=(K0) + |0,n)): e )>(-L-J) 2 . 
V2 n 7r 

The lower bound (— [^rj ) 2 converges e 2 . This means that 
the mini-max of the mean square error does not decrease 
when n increase. Hence, 

lim n 2 Cg (-^(\n,0) + \0,n)):e)=oo. (20) 
n-»oo y/2 

Next, we consider the asymptotic limit of the per- 
formance with the global asymptotic min-max criterion 
when the noon state i(|n, 0) + |0, n)). By using (|20|) 

and ((7]), the relation 

lim n 2 C(-^(\n,0) + |0,n») = oo (21) 

holds. This fact can be shown in another method based 
on (|13[) . In fact, the L 2 function corresponding to the 
noon state via (|T0|) is fo(x) := ^{S(x — 1) + S(x + 1)). 

The function f n (x) := -±{5{x ~ 1) + e"> 5(x + 1)) satisfies 

(fn\ p2 \fn) = oo- Thus 5 the relation (PJ yields (j2T|) . 

Next, we replace the number state by the coherent 
state. In the following, the coherent state with complex 
amplitude «i,a2 is written as |ai, 0:2)0 Then, we con- 
sider the case when the covariant measurement is applied 



and we input the noon state 4j(|n, 0) c + |0,n) c ), which 
is called the coherent noon state. The L 2 function corre- 
sponding to the coherent noon state via (fTU|) is the same 
function /q. Hence, the relation (fT3]) yields 



lim n 2 C(^(\n, 0) c + |0, n) c )) = 00. 

n->oo y/2 



(22) 



VII. CONCLUSION 

We have discussed three kinds of bounds of MSE in 
the phase estimation with photon number constraints. 
The first is the global minimax bound. The second is 
the local minimax bound. The third is the bound with 
the locally unbiased condition. These bounds have been 
treated with two kinds of constraints: One is the con- 
straint for the average of the square of the photon number 
concerning the input state. The other is the constraint 
for the maximum photon number concerning the input 
state. We have shown that the asymptotic limits C av , 
C aVt $, and C av ,e of three kinds of bounds coincide un- 
der the first constraint. However, these bounds C max , 
Cmaxfii and Cmax, 9 do not coincide under the second 
kind of constraint. In fact, the locally unbiased condition 
and Fisher information are originally mathematical con- 
cepts. In the independent and identical distributed set- 
ting of state estimation, these values provide the asymp- 
totic bound of MSEpjj-Qjl- However, these values do 
not provide the operational meaning in general. That is, 
these values have no operational meaning when they do 
not coincide with operational values. Hence, under the 
second kind of constraint, only these bounds C max and 
C m ax,e have the real meaning. 

Further, we should be careful of the meaning of the 
noon state and the coherent noon state. When these 
states are inputcd, the optimal SLD Fisher information 
is realized. However, under the local mini-max criterion, 
as is discussed in section 6, the MSE of the noon input 
state does not convergence to zero, and the noon input 
state is far from the optimal input. The optimal input is 
given in and (fT5)l . Hence, in order to realize the high 
performance phase estimation, it is desired to implement 
the input states and (fT5|) . 
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